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Analyzing recent experimental results 1,2 , we find similar behaviors and a deep analogy between 
three-junction superconducting qubits and suspended carbon nanotubes. When these different sys- 
tems are ac-driven near their resonances, the resonance single-peak, observed at weak driving, splits 
into two sub-peaks (Fig. 1) when the driving increases. This unusual behavior can be explained 
by considering quantum tunneling in a double well potential for both systems. Inspired by these 
experiments, we propose a mechanical qubit based on buckling nanobars — a NEMS so small as to 
be quantum coherent. To establish buckling nanobars as legitimate candidates for qubits, we cal- 
culate the effective buckling potential that produces the two-level system and identify the tunnel 
coupling between the two local states. We propose different designs of nanomechanical qubits and 
describe how they can be manipulated. Also, we outline possible decoherence channels and detection 
schemes. A comparison between nanobars and well studied superconducting qubits suggests several 
future experiments on quantum electromechanics. 

PACS numbers: 85.85.+j 



I. INTRODUCTION 

Micro- and nano-electromechanical systems (MEMS 
and NEMS) have attracted widespread attention because 
of their broad spectrum of functionalities, their tiny sub- 
micron sizes, and their unique position bridging micro- 
electronic and mechanical functions 3 . Sophisticated tools 
ranging from mirrors and sensors, to motors and multi- 
functional devices have been fabricated 4,5 . As the size 
of these devices shrinks, experimental studies of NEMs 
are rapidly approaching the quantum limit of mechanical 
oscillation a 5 i 6 i 7 i 8 i 9 , where quantum coherence and super- 
position should result in quantum parallelism and the 
possibility of information processing. 

The emergence of quantum electromechanical devices 
(see, e.g., Jief&^&Z&l and references therein) brings 
both challenges, such as the inevitable and ubiquitous 
quantum noise, and promises, such as macroscopic quan- 
tum coherence^iiSiii or quantum teleportationi 2 * 1 ^ . In- 
deed, during the past several years the quantum mechan- 
ical properties of NEMS and how they can be coupled to 
other quantum mechanical objects have been very ac- 
tively studiedi^ii^iiSiiiiSiiS. The experimental pursuit 
of these studies has so far focused on cooling a device to 
reach states where quantum fluctuations in the lowest en- 
ergy eigenmode dominate over thermal fluctuations^^ 5 -. 
Such eigenmodes are generally harmonic-oscillator modes 
with equal energy-spacings and follow bosonic statistics. 

Among NEMS there also exist systems that can be well 
approximated by two levels in the quantum limit, so that 
they might be candidates for qubits. For example, when 
a longitudinal strain above a certain critical value is ap- 



plied to a small bar with one or two ends fixed, there exist 
two degenerate buckling modes, as schematically shown 
in Fig. 2. In the quantum mechanical limit, these two 
modes represent the two lowest-energy states. They can 
be well separated from the higher-energy excited states, 
so that at low temperatures the buckling nanobar can be 
properly described as a two-level system. Furthermore, 
since nanobars can be charged or can carry electrical cur- 
rent, the electric or magnetic field can be conveniently 
used to manipulate their quantum states. 



The fascinating prospect of observing quantum coher- 
ent phenomena in a macroscopic mechanical oscillator is 
a main motivation of this study. To achieve this goal, a 
NEMS needs to possess at least two general attributes: 
small size and a high fundamental frequency. A carbon 
nanotube is a natural candidate: it is very thin (and can 
be very short) while still being stiff 20 (thus having high 
vibration frequency) because of the strong covalent bonds 
between carbon atoms within the graphene sheet^i*2£. 
However, without strong experimental support, it might 
be hard to judge the feasibility of achieving macroscopic 
quantum coherence in NEMS. In contrast to NEMS, ex- 
perimental evidence has demonstrated that supercon- 
ducting charge^ 3 -, flu^4, and phased qubits, using Su- 
perconducting Quantum Interference Device (SQUID), 
do exhibit quantum coherence. Thus, it is important to 
look for possible similarities in the behavior of super- 
conducting qubits and NEMS. This should help future 
prospects to observe quantum coherence in NEMS. 



2 



II. SPLITTING RESONANCE AS A 
MANIFESTATION OF QUANTUM 
TUNNELLING IN THREE-JUNCTION SQUIDS 
AND NEMS 

Experiments 1 on suspended single-wall carbon nan- 
otubes (with diameter of about 20 nm and length of 
about 1.7 /im) excited by an electromagnetic wave show 
a resonance peak with an unusual shape (Fig. la), for 
one of the fundamental harmonics uoq. As expected, a 
Lorentz-form resonance peak, that grows with increas- 
ing intensity of the ac field, was observed at weak driv- 
ings. Surprisingly, when further increasing the ampli- 
tude of the externally applied electromagnetic wave, this 
peak splits into two sub-peaks 1 . When further increas- 
ing the ac drive, these two sub-peaks gradually move 
away from each other, while their heights stop growing. 
It is important to stress that this phenomenon was ob- 
served at a frequency (o;o/27r ~ 2 GHz) and temperature 
(T ~ 100 mK) where quantum effects start to domi- 
nate over thermal noise (T < hcoo/kB ~ 100 mK, with 
Boltzmann constant Also the dissipation in the sys- 
tem was quite low (quality factor Q ~ 1500), which 
is important to observe quantum effects. Note that the 
fundamental frequency of this device can be easily in- 
creased (at least by an order of magnitude, well into the 
operating frequency 1-15 GHz of many superconducting 
qubits), for experimentally available carbon nanotubes 
with shorter length. 

Interestingly, a similar phenomenon has been recently 
found 2 for an Al three-junction SQUID qubit coupled 
to a Nb resonant tank circuit (Fig. lb). It was exper- 
imentally proven 2 , via the observation of quantum hys- 
teresis (Landau- Zener transitions), that this circuit was 
operated in the quantum regime (at coo/2tt ~ 20 MHz 
and T > 10 mK), though for a worse ratio of quan- 
tum to thermal noise compared to the carbon nanotubes^ 
(i.e., coo/T ~ one order of magnitude higher for the nan- 
otube). When the magnetic flux in the SQUID was 
driven as $ = <l>dc + $ac cos(cjt), the resonance in the 
response, probed via the tank voltage as a function 
of the dc flux, was found 2 to exhibit a transformation 
from a Lorentz-form single-peak to a double-peak shape 
(Fig. lb) in striking similarity to suspended driven car- 
bon nanotubes 1 . Sweeping the dc flux in the SQUID 
corresponds to changing the fundamental frequency as 

ujo = ^o($dc = 0) + Acj ($dc), (1) 

with Aujq oc (<E>dc — $o/2) and flux quantum <3>o- Thus, 
we find that the measured 1,2 response, near resonance, 
of both systems (carbon nanotube and three-junction 
SQUID qubit) depends on the difference uo — uoq. This 
will be useful to establish below that the measurements 
for driven carbon nanotubes^ and three-junction SQUID 
qubit 2 essentially probe the same effect. 
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FIG. 1: (a) Experimental resonance- of the effective dissi- 
pation (which is proportional to the average energy (E)) as 
a function of the dimensionless frequency (u — cjo)/<^o, for 
a suspended carbon nanotube, for different intensities P of 
the externally applied ac electromagnetic wave, (b): the res- 
onance in the response of the three-junction SQUID probed 
by measuring 2 the voltage for the resonant tank Vr for dif- 
ferent driving amplitudes as a function of (oo — ooo) / ooo . For 
these different systems, a similar splitting of the resonance 
peak was observed^. Schematic diagrams for a driven sus- 
pended carbon nanotube and three junction SQUID are also 
shown at the top. The simulated resonance [average energy 
(E) versus the dimensionless frequency {uo — ujq) / uoq] for quan- 
tum (c) and classical (d) particles moving in the double-well 
potential shown in the inset of (c). A quantum description 
agrees with experimentsi^Fig. la,b): For quantum particles, 
the standard Lorentz-form resonance peak for weak ac drives 
(c, red curve, A = 0.005) splits into two sub-peaks due to 
tunnelling at stronger drives (blue curve, A = 0.01; green 
curve, A = 0.0135). A classical description cannot describe 
experiments 1,2 (Fig. la,b): For classical particles, the weak- 
driving Lorentz-form peak (red curve, A = 0.005) becomes 
asymmetric and exhibit discontinuous jumps E(oo) (green 
curve, A = 0.0135) due to the nonlinearity of U(y). 
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We propose how the incorporation of quantum tun- 
nelling can be used to understand these still-unexplained 
experimental observations in nanotubesA and three junc- 
tion SQUIDs 2 . In order to interpret the splitting of the 
resonance peak (Fig. la and lb), we employ a quasi- 
classical model 26 described by the equation of motion 



y + 2Xy + dU(y)/dy = A cos(c^), 



(2) 



for the double- well potential shown in the inset in Fig. lc, 
with small damping A = 0.01, and different driving am- 
plitudes A. We assume that the particle (which mimics 
the buckling mode or current in the SQUID) can tunnel 
from one well to the other with transmission coefficient 



D = exp 



2V2 



/ VTT^Edy 

J —a 



(3) 



where a is the classical turning point, and E = (y) 2 /2 + 
U(y) is the energy of the particle. This minimal model 
allows to qualitatively describe the resonance of a quan- 
tum particle in a double well potential (a more com- 
plete, fully quantum mechanical, theory will be published 
elsewhere 27 ) and also to obtain the transition to a clas- 
sical description (Fig. lc,d). 

When the potential barrier Uq is comparable with fiwo, 
we find that the single resonance peak, at low driving, 
splits into two sub-peaks (Fig. lc) for higher drives; in 
agreement with experimental findings 1,2 . The physical 
origin of this effect is the following: (i) At low driv- 
ing the energy E of the quantum particle (whether the 
SQUID or the nanotube) is also low and the probability 
of tunnelling is negligibly small, thus, the usual Lorentz- 
form of the resonance occurs; (ii) When the driving (and, 
thus, energy) increases, the particle starts to tunnel be- 
tween the wells. A tunnelling event is equivalent to an 
inversion around the origin, i.e., changing y to —y: like 
replacing the left well by the right one, or vice versa. 
As a result, right after a tunnelling event, the average 
power V — (yF cos cot) the particle receives from the ex- 
ternal force (e.g., the power absorption by the nanotube 
from the external electromagnetic wave) changes its sign 
V — > — V and, hence, the energy, which increases be- 
fore tunnelling, starts to decrease temporarily. Thus, the 
average energy (E) also decreases when the number of 
tunnelling events increases. This explains the minimum 
of (E(co)) at uj = coo (instead of the maximum observed at 
weak driving) and also the splitting of the resonance peak 
- because the tunnelling frequency sharply increases at 
the resonance frequency uoq. 

For the case of "more classical" particles, Uq ^> Tiujq, 
the probability of tunnelling is always very low. When 
driving increases, the particle (which is always located in 
the same well) begins to feel the strong nonlinearity of the 
potential U (?/), resulting in the energy dependence of the 
oscillation frequency ljq(E). Instead of a split resonance, 
the resonance peak as a function of frequency shows the 
standard asymmetric shape with a sharp jump of (E(uj)), 
associated with mechanical hysteresis 28 (Fig. Id). 



Therefore, we conclude that the two-peak resonance in- 
dicates that both the three-junction SQUID qubit 2 and 
the suspended carbon nanotubeA, operate in the quan- 
tum regime. Thus, the technology for fabricating sus- 
pended buckled carbon nanotubes, working as nanome- 
chanical qubit s, already exists 1 . In view of the explo- 
sive growth of NEMS technology, below we discuss the 
prospect of such buckling charged nanobars (the clamp- 
ing at the base ensures an anisotropic nanobar instead 
of an isotropic nanotube) as candidates of quantum bits 
for quantum information processing. Indeed, previous 
work on quantum buckling (e.g.^Ai) were not focused 
on operating these as qubits, which is our focus here 2 ^. 



III. CONTROLLING QUANTUM STATES OF A 
NANOROD IN A DOUBLE- WELL POTENTIAL. 



The double-well potential (Fig. 2a) corresponding to 
the buckling modes y sin[7r/// max ] of a nanorod (Fig. 2b) 
takes the form: 



U(y)=ay 2 +Py 4 + 



2lm&xf± 



-y, 



(4) 



where the parameters can be estimated^ as 

^2 



<*(/) = 

/3(f) = 

fc = 



Al 



ifc-f), 
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4 



64^ max 
IYtt 2 
I 2 

max 



(4/c-3/), 



(5) 



Here, we introduce the tube length I (0 < I < /max), the 
elastic modulus Y and the moment of inertia /. To con- 
trol the rod we need: (i) a longitudinal compressing force 
/ acting on the rod ends, and (ii) a transverse force f±, 
which can be produced via, e.g., interacting the charged 
nanorod with an electric field. 

For zero transverse force f± = 0, our qubit (Fig. 2b) 
is in its degeneracy point. Indee d, thi s potential has two 
minima, at y = ± 2/0 (/) = ±vW^ that are separated 
by a potential barrier 



AU(f) = a 2 /4p. 



(6) 



The first two energy levels E\ and E2 in the right well can 
be estimated assuming a parabolic potential well shape 
U w mujQ^y — yo) 2 /2 with 



m 



m 



1/2 



(7) 



and the mass m of the nanorod. Thus, we obtain 

E 1 (/) = 3a(a(/)/m) 1 / 2 (8) 

and 

E 2 (f) = 5E 1 (f)/3. (9) 



Due to quantum tunnelling between the left and right 
potential wells, E\ splits into two levels 



Ef = E x ± KA t . 



(10) 



The tunnelling rate At between the left and right buckled 



states is: 



Mf) 



.31 



«(/) 



■ exp 



/2^(AU(f) - E.jf)) 



(11) 

It is important to stress that the longitudinal force / 
allows to control the tunnelling A t (f) as well as the en- 
ergy levels itself E\(f). The higher levels are well sepa- 
rated from the two lowest ones: (E± — E^)/(E2 — E\) ~ 
At/hujQ <C 1. Changing the transverse force f± moves 
the system out of the degeneracy point, allowing to ma- 
nipulate the proposed nano- mechanical qubit. 

We now introduce another, more controllable, design 
of a nanorod qubit, which is also less affected by deco- 
herence coming from the relative vibrations of the rod 
ends (Fig. 2c). In such a design we consider a "crossed" 
electric field having a component E\\ along the rod and 
another E± perpendicular to the rod, which is clamped 
to a substrate only by one end. Changing the charge of 
the top capacitor plate and the charge of the rod itself 
can control both the electric field, E\\, on the rod and 
its gradient, dEu/dx. By adjusting the different knobs, 
we can tune the relative strength between thermal acti- 
vation and quantum tunneling, allowing the observation 
of transition between these two regimes. 



IV. MANIPULATION, DECOHERENCE, AND 
DETECTION OF MECHANICAL QUBITS 

The manipulation of the mechanical qubits can be 
achieved electrically (Fig. 3). For example, in analogy 
to the Cooper pair box 23 (see table), one can prepare the 
nanobar qubit in the \L) state by setting a transverse 
electric field towards the right (assuming the nanobar is 
negatively charged). By suddenly turning off this electric 
field and bringing the system to the degeneracy point, the 
nanobar state is prepared in a coherent superposition of 
{\R) ± \L))/y/2. Because of the tunnel splitting, the sys- 
tem then starts to oscillate coherently, with a frequency 
given by A t . Therefore, by detecting the nanobar posi- 
tion, as a function of At (see Fig. 3), one can determine 
the coherent oscillation frequency and the system deco- 
herence. Driven coherent transitions between the two 
qubit states can also be similarly achieved. A sinusoidal 
component can be added to the vertical electric field that 
is used to control the tunnel coupling between the \L) 
and \R) states, in analogy to the microwave driving force 
on the Josephson coupling in a Cooper pair box. The 
study of free and driven coherent oscillations of a nano- 
bar would help clarify its quantum coherence properties 
and demonstrate its feasibility as a qubit. 





(b) 




FIG. 2: (a) Double-well potential for a buckled nanobar. Due 
to tunnelling from the right potential well to the left one, 
the lowest energy level is split into two levels for f± = as 
shown in the left panel. The lowest (blue) and the first excited 
(red) levels correspond to the symmetric and antisymmetric 
combinations of the wave functions localized in the left and 
right potential wells. The energy level splitting between the 
left and right states can be controlled by the transverse force 
f± as shown in the top right panel, (b) A buckled rod qubit 
where the compressed force applied to the rod ends controls 
the potential shape [a and (3 in Eq. and, therefore, the 
energy splitting at the degeneracy point. The transverse force 
f± allows to drive the bar to a degeneracy point, (c) Another 
proposed design for a buckled-rod qubit. In this case the 
possible decoherence originating from the relative vibration of 
the top and bottom rod holders can be avoided. By changing 
the charge of the top capacitor plate and the charge of the rod 
itself it is possible to independently control both the parallel- 
to-rod electric field E\\ and its gradient. This design allows 
a high level of control of both the tunnelling and the energy 
splitting, (d) Correlated noise produced by some phonons can 
be eliminated using a decoherence- free subspace: qubits Ql 
and Q2 can be associated with one logic qubit. (e) Single- 
shot measurements could be done via either single electron 
transistor (SET) or photon reflection. 



Universal quantum computing requires two-qubit op- 
erations. For charged nanobars, the inter-qubit interac- 
tion comes naturally in terms of the electric dipole in- 
teraction between the bars, quite similar to the dipole 
interaction used in other proposed qubits. In the case 




time 



FIG. 3: A schematic diagram for controlling the quantum 
state of a nanobar via coherent oscillations at the degener- 
acy point. Tuning the system to its degeneracy point can 
be done by changing the perpendicular electric field E±. For 
the case when E± and E\\ is zero at t = 0, the perpendicular 
electrical field E± has to be applied first and then the longitu- 
dinal field E\\ (e.g., at the moment marked by the solid circle 
on the E±(time) curve). This initializes the system during 
the stage (i). Switching off the perpendicular electric field 
E± brings the system back to its degeneracy point, and the 
nanorod starts to oscillate a time At during stage (ii). To ob- 
serve coherent oscillations, measurements (e.g., optically or 
electrically) must be done for many values of At. 



of nanobar qubits, one can again use the transverse elec- 
tric field Et to tune selected qubits into resonance, then 
apply microwaves to perform conditional rotations and 
other operations. 

For NEMS, and for more conventional applications 
such as resonators, major sources of noiso 32 ! 33 (that limit 
the quality factor Q) are internal t her mo mechanical noise 
(such as heat flows due to inhomogeneous distribution of 
strain), Nyquist- Johnson noise from the driving electrical 
circuit, adsorption-desorption noise when a resonantor is 
moving in a non- vacuum, noise from moving defects, etc. 
In the quantum mechanical limit, some of these noises be- 
come unimportant. For example, there should be no heat 
flow: the nanobar should not be driven beyond its first 
excited state. The NEMS can be placed in a vacuum to 
reduce adsorption-desorption noise. The whole system, 
including its electrical components, has to be cooled in 
a dilution refrigerator in order to reach the two lowest 



System 


SC charge QB 


SC flux QB 


Nano-bar QB 


States 


I excess charge > : 

| ) and | 1 ) 


1 current direction > : 

IC> md 15) 


1 buckling direction > : 

| R ) and | L ) 


Hamiltonian 


H = e a z + h A t cr x 


Tunneling £±t and 
energy splitting 
controlled by 


Gate voltage 

(normalized: Tig ) 


Magnetic flux: c[) 

(normalized: 

/ = $/$o ) 


Transverse force /j_ 

(induced, e.g., by a 
transverse electric 
field E± ) 


Tunneling 
controlled by 


Josephson energy: 


Josephson energy: 

Sj 


Either 

• Longitudinal force 

(induced by pressure or 
electrical field fu ) 

• or its gradient 
Longitudinal force and 
its gradient control d x /i 
barrier height and its 
curvature 


Coupling 
between qubits 


Electrical 

(e.g., capacitive or 
inductive coupling) 


Magnetic 


Electrical (dipolar) 


Decoherence 
sources include 


Charge fluctuations 


Flux fluctuations 


1. Charge fluctuations 

2. Phonon-phonon 
interactions 


Read-out 


Electrical 

(e.g., SET or JJ) 


Magnetic (SQUID) 


Either electric or 
optical, or mechanical 



Table: comparison of Josephson-j unction superconducting 
(JJ SC) charge, JJ SC flux and nano-bar qubits (QB). 

states of the nanobar, so that the Nyquist- Johnson noise 
is also suppressed. The main source of quantum me- 
chanical decoherence might be internal dissipation caused 
by phonon-phonon interactions. Since the nanobars are 
clamped onto a much larger substrate, in which lower 
energy phonon modes exist, coupling at the base of the 
nanobars could lead to relaxation/excitation of the qubit 
states. However, the long- wave- length phonons produce 
a correlated noise, which can be reduced by working 
within a properly chosen decoherence-free subspace (see 
schematics in Fig. 2d and, e.g.j2L2A). Another source 
of decoherence would be charge fluctuations. Charged 
nanobars allow increased maneuverability of the NEMS; 
but the surrounding and internal charge fluctuates can 
produce charge noises that affect the relaxation and de- 
phasing of the qubit. 

Single-wall nanotubes are a natural candidate for the 
mechanical qubit we consider here. However, our pro- 
posal is not limited to them. Indeed, a multiple- wall nan- 
otube could provide a higher-frequency mechanical oscil- 
lator and more favorable condition for observing macro- 
scopic quantum tunnelling and the coherent evolution of 
mechanical motion. Silicon-based systems provide an- 
other enticing alternative, especially from the perspective 
of fabrication. While it might be difficult to clamp several 
single-wall nanotubes to make identical nanobars (with 
controlled inter-bar distance, same length, same buckling 
orientation, and same response to external stress), fabri- 
cating a lithographically-patterned silicon-based nanos- 
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tructure would be much more reliable. SiC nanobars 
have shown higher stiffness 22 compared to Si. Unless 
nano-assembled nanotubes can be made with sufficiently 
high precision, materials that can be lithographically fab- 
ricated seem more promising candidates for a larger-scale 
mechanical quantum information processor. This would 
be an ironic turn of events, given that the first computers 
(by C. Babbage) were mechanical. 

The detection of the nanobar state can be done either 
electrically2ii2iS or optically 36,37 . It has been shown in 
recent experiments^ that single electron transistors are 
very sensitive to small charge displacements. Optical de- 
tectors can also be used, where light scattering can detect 
the state of the bent nanorod (Fig. 2e). 

Very recently, Ref. 38 fabricated suspended nanobars 
(driven by a 25 MHz current through an attached elec- 
trode) switching between two distinct states. These sus- 



pended nanobars have already been tested 3 ^ as very fast 
and very low-power-consumption storage memory de- 
vices. Still, many challenges lie ahead on the road to 
practical quantum electromechanics. We hope that our 
proposal here stimulate more research towards the ulti- 
mate quantum limit of NEMS. 
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